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ABSTRACT 

We perform detailed magnetohydrodynamic simulations of the gas cloud G2 interacting with 
the accretion flow around the Galactic Center black hole Sgr A*. We take as our initial con- 
ditions a steady-state, converged solution of the accretion flow obtained earlier using the 
general-relativistic magnetohydrodynamic code HARM. Using the observed parameters for 
the cloud's orbit, we compute the interaction of the cloud with the ambient gas and identify 
the shock structure that forms ahead of the cloud. We show that for many configurations, the 
cloud front crosses orbit pericenter 7 to 9 months earlier than the center-of-mass. 

Key words: accretion, accretion disks, black hole physics, relativity, acceleration of particles, 
radiation mechanisms: non-thermal 



1 INTRODUCTION 



iGillessen et alJ J2012L 1201 3l) discovered an object, likely a cloud 
of molecular gas, called G2 that is approaching the Galactic Cen- 
ter on a highly eccentric orbit. The center-of-mass of the cloud is 
expected to reach pericenter, located at radius R « 3.0 x 10 15 cm 
from the supermassive black hole Sagittarius A* (Sgr A*), at to = 
2013.69 (September 2013fl During this encounter, the cloud will 
move supersonically thr ough the hot accretion flow around Sgr A* 



Sadows ki et al. 12013b) and a bow shock will form ahead of it 
Naravan et alJ I l2012ah . The shock provides suitable conditions for 
particle acceleration and is likely to accelerate hot ambient elec- 
trons itordatMstic energies . 

iNaravan et all d2012ah and lSadowski et all d2013bh estimated 
the efficiency of this acceleration and calculated the expected ra- 
dio synchrotron emission. Depending on model parameters, they 
predicted radio fluxes ~ 1 - 20 Jy at v = 1.4 GHz. In addition, us- 
ing simplified assumptions, Sado wski et aDd2Q13bh estimated that 
the bow shock will cross pericenter roughly 5 months before the 
center-of-mass does, and that the peak radio flux will be reached 
roughly a month later. However, their estimate did not take into 
account the structure of the cloud or the details of the dynamical 
interactions between the cloud and the ambient gas, both of which 
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1 Recently, Phifer et alJ d20 1 31) obtained another orbital solution, derived 
from Br-gamma line astrometry and radial velocity measurements, resulting 
in a later time of closest approach (to = 2014.21, March 2014), slightly 
closer periastron (3800 vs 4400Rq), and a longer orbital per iod (276 vs 198 
yrs). I n this work we use the orbital parameters derived by IGillessen et al.l 
(2013). Because of the similarity of the orbits the shock parameters and its 
location with respect to the center-of-mass are expected to be similar. 



have the potential to change the shock properties and the expected 
radio signals. 

In this paper we obtain a more reliable estimate of the shock 
location and the interaction epoch by carrying out general relativis- 
tic magnetohydrodynamic (GRMHD) simulations of G2 moving 
through the turbulent hot accretion flow around Sgr A*. The Simula 



tions are performed with the GRMHD code HARM (Gammie et al 



2003) using a realistic model of the accretion flow (Narava n et al 



l2012bl) .The cloud is scaled down to fit within the converged region 
of the disk simulation. Although the spatial resolution of the simu- 
lations is limited, we are able to determine with reasonable confi- 
dence that the bow shock is formed ahead of the cloud, and that it 
crosses pericenter 7 to 9 months earlier than t . The accompanying 
radio emission is thus exp ected to reach peak va lues in Spring or 
late Su mmer 2013 for the IGillessen et all d2013h and IPhifer et al.l 
b013h orbital models, respectively. 

Our simulations are the first that include magnetic fields and 
propagate the cloud through a turbulent accretion flow in three di- 
mensions (3D), thus probing phenom ena that were not explored in 
previous studies. iBurkert et al.l d201 2l) and lSchartmann et al.ld2012h 
performed two-dimensional fixed-mesh simulations of the cloud 
approaching Sgr A* and penetrating an isotropic, smooth, non- 
magnetized and artificially stabilized atmosphere. Si milar assump- 
tions about the background gas were adopted by lAnninos et al.l 
d2012h who performed 3D adaptive mesh simulation, which are 
best suited for resolving instabilities at the cloud surface (however, 
they s till neglected the role of magnetic fields). Finally. [Saitoh et all 
(2012) used a Smoothed Particle Hydrodynamics (SPH) code to 
simulate the evolution of the cloud assuming the same model for 
the ambient gas as the other authors, but allowing the gas to cool 
through optically-thin radiation. 

The structure of the paper is as follows. In Section we dis- 
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cuss the initial conditions for our simulations, in Section [3] we de- 
scribe the simulation results, and in Section [4] we discuss uncer- 
tainties and the optimal strategy for radio observations of the G2 
impact. 



2 MODEL 

2.1 Disk model 

The radiatively inefficient accretion flow around Sgr A* has been 
studied extensively through observations and theoretical investiga- 
tions. The best available constraints come at large radii from the 
measurement of the gas density at the Bondi radius, which yields 
n e (Rz) ~ 130 cm -3 (Baganoff et al. 2003), and in the innermost re- 
gions, where the mass accretion rate onto the bl ack hole has been 
estimated to be M RH ~ 10~ 9 + 10~ 8 M yr" 1 (see lYuan et al|2003| : 
Marrone eta"il2 007: Moscibrodzka et al. 20091: iDexter et alJboiol : 
Shcherbakov et al.ll2012l : iDibi et al.ll2012h . Models of the gas den- 
sity and the electron temperature at these small distances can simul- 
taneously account for the long-wavelength spectrum of the accre- 
tion flow, the polarization measurements, and the size of the emit- 
ting region at 1.3 mm (see Broderick & Loeb 2005). Less is known 
directly about the properties of the flow at intermediate radii, where 
the pericenter of G2's orbit lies. Psaltis (2012) carried out a compar- 
ison of the gas properties obtained at intermediate radii by extrapo- 
lating the flow properties inwards from the Bondi radius using the 
model of Quataert (2004) and outwards from the model of Broder- 
ick & Loeb (2005) and found reasonable agreement between these 
solutions between few xlO 2 to ~ 10 5 /?g where Rq = GM/c 2 and 
M is the mass of the black hole. 

Numerical simulations similarly rely on the measurements at 
small and large radii as anchor points to calculate the properties 
of the gas throughout the flow. In this paper, we make use of the 
radiatively inefficient a ccretion simulation aro und a non- spinning 
black hole described in Narava n et al . (2012b), which corresponds 
to the so-called Magnetically Arrested Disk (MAD, Narayan et al. 
2003) regime. The accretion disk was solved on a fixed, horizon- 
penetrating grid with 264 cells logarithmicaly spaced in radius, 126 
cells roughly uniform in 6 and 60 cells uniform in 0. It was initial- 
ized as an equilbrium to rus threaded by seed, single-loop magnetic 
field JPenna et al.ll2012h . Magneto-Rotational Instability (MRI) de- 
velops, the gas accretes and it brings significant magnetic flux onto 
the BH. We adopted an ideal equation of state with adiabatic index 
y = 5/3. Fig. Q] shows a typical snapshot of the gas density in the 
r6 plane. 

Because of the very long duration of this simulation, t = 

2.2 x 10 5 R G /c, the simulated accretion flow is in steady state out 
to a relatively large radius ~ 200R G at the equatorial plane, and to 
much larger radii at higher latitudes. We normalize the gas density 
in the simulation such that it is consistent with the measured den- 
sity at the Bondi radius and, at the same time quantitatively agrees 
with the expected accretion rate at the BH horizon. Shcherbakov 
(2013, private communication) have shown that the simulated spec- 
trum based on this model agrees with observed spectrum of Sgr A*. 
Fitting fo rmulae for the mean properties of the accretion flow are 
given m Sadowski etal. (2013b). The mass of the disk gas inside 
the cloud radius of pericenter is roughly 0.1M Ear th and is much 
small er than the mass of the cloud M c \ « 3M Earth dGillessen et al.l 
l2012h . 




x(R G ) 



Figure 1. Snapshot of density (arbitrary units) in the r6 plane for the back- 
ground model (evolved without the cloud) of the accretion disk taken at the 
same time as in Fig. [3] (which shows the results of cloud and disk interac- 
tion). 

2.2 Cloud model 

To set up the configuration of the cloud in phase space, which we 
use as initial conditions for our cloud- disk interaction s imula tions, 
we follow the approach described in iGillessen et al.l (12012b . We 
start with a spherical cloud of test particles at an earlier epoch cor- 
responding to year 2000.0. We adjust the location and the velocity 
of the center-of-mass of th is initial cloud so that they are consistent 
with the orbit described in IGillessen et al.l (120131) . We also assign a 
spherically symmetric density profile and a distribution of random 
velocities to the cloud (described in detail in the next paragraph) 
and adjust these so as to satisfy the distributions o f positions and 
velocities observed in year 2012 (Gillessen et al .120131) . This proce- 
dure does not yield a unique model f or the initial cloud. Th erefore, 
in addition to the model described in Gillessen et al . (2012), which 
we refer to as model Nl, we consider a number of other models 
which are all consistent with the data (see Table [TJ- Our aim is to 
explore a range of models in order to assess the effect of the ini- 
tial cloud configuration on the location and properties of the bow 
shock. 

In the alternative models we consider, we model the density 
profile and the velocity distribution using two parametrizations for 
the former and three for the latter. For the initial density profile 
(year 2000), we consider a Gaussian distribution, with two values 
of the full width at half maximum (FWHM), and a constant density 
cloud (called flat), with three choices of the outer radius R out (Ta- 
ble [J}. For the Gaussian density model, we truncate the cloud at a 
radius equal to 2<x, where cr is for the standard deviation. For the ve- 
locity distribution, we take: 1. a Gaussia n distribution with FWHM 
= 120 km s _1 as in lGillessen et al.ll2012l (called normal), 2. a distri- 
bution where the magnitude of the additional velocity component is 
proportional to the distance from the center-of-mass, with a veloc- 
ity of 120 kms -1 at FWHM or R out for Gaussian and flat models, 
respectively (called scaled), and 3. a model with no velocity dis- 
persion, where all the particles in the cloud move with the same 
velocity as the center-of-mass (called no). We give in Table Q] the 
details corresponding to all the models we have considered. Each 
of these models, when propagated forward from year 2000 to year 
2012, gives a cloud of particles with positions and velocities that 
are consistent with observations. 

2.3 Initial MHD setup 

The GRMHD s imulation that we u se to model the background ac- 
cretion flow dNaravan et al.ll2012bl) has achieved steady state out to 
a radius R ~ 200R G at the equatorial plane, whereas the cloud orbit 
has a pericentric radius of 4400/? G - To evolve the cloud-disk system 
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numerically we need to either extrapolate the disk model to several 
thousand R G or scale down the cloud and its orbit to ~ IOO^q- We 
choose to do the latter since it allows us to retain the full 3D turbu- 
lent structure of the simulated disk. 

We scale all distances down by a factor X = 40, which brings 
the pericenter of the orbit down to R = HO^g- This choice, al- 
though "ad hoc", ensures that, on the one hand, the cloud fits within 
the converged region of the simulation and, on the other hand, it 
does not go too close t o the black hole, where the disk model is no 
longer self- similar fsee lSadowski et al.ll2013bf) . We scale the cloud 
size, orbit semi-major axis, and other relevant distances by a factor 
of X~ l , the orbital velocity by X 1/2 , and the time relative to pericen- 
ter by X~ 3/2 . The gas densi ty in the simulation varies with radius as 
R~ l dSadowski et al . 2013b). Therefore, we scale the cloud density 
also by X" 1 . 

The strength of the expected bow shock depends on its Mach 
number — the ratio of the relative cloud-disk velocity to the disk 
sound speed. The relative velocity increases proportional to the 
orbital velocities, i.e., by X 1/2 . The temperature in the numerical 
model follows X 1 and therefore the speed of sound scales as X 1/2 . 
As a result, the shock strength is independent of the adopted value 
of the scaling factor X. 

We report all simulation results in length units of R G = 
XGM Sgr A*/c 2 , though when reporting times, e.g., in the Tables, we 
convert back to the parameters of the original G2 cloud. 

As the initial conditions for the disk, we choose the final con- 
figuration ayimej_=2.2j<J0 5 ^G/c of the ADAF/MAD run de- 
scribed in lNaravan et al.ld2012bh . Then, for each of the cloud mod- 
els shown in Table [TJ we evolve the test particles forward along 
geodesies up to time t mM = 2012.5 (except for one model, NOF, 
where we choose time = 2012.0; see below). We transfer the 
mass and momentum corresponding to each test particle to the cor- 
responding grid cell in the simulation, thereby adding the cloud 
density and momentum to the pre-existing accretion gas in that cell. 
We assume that the cloud gas is cold (this is certainly true relative 
to the > 10 8 K ambient gas in the accretion disk) and unmagne- 
tized. Thus, the pressure and magnetic field retain the same val- 
ues as in the original GRMHD simulation. Starting from this initial 
state (Fig. [3, we run the simulation forward in time and study the 
motion of the cloud as well as its interactions with the medium. 

We ran additional models to test the effects of the initial setup 
on our results. We ran model NOF with a start time of 2012.0 and 
confirmed that selecting a start time of 2012.5 vs. this earlier time 
does not modify the results. In order to assess the impact of the 
ambient magnetic field on the cloud motion, we also ran one model 
(NOB) which is identical to the fiducial model NO except for a mag- 
netic field strength that we set equal to zero. We will discuss the 
results of this latter run in the later sections. 

2.4 Orbit orientation 

The GRMHD simulation of the background accretion disk covers 
the whole domain of (r, 0, 0) and has sufficient resolution (264 grid 
cells in r, 126 in 0, 60 in 0) to capture the physics of the domi- 
nant MRI modes. The cloud, however, has a relatively small ex- 
tent compared to the full domain. Therefore, the simulations we 
describe here have only modest resolution on the scale of the cloud; 
in particular, our resolution is coarser than other simulations of G2 
descr i bed in the literatu re dAnninos et al . 2012; S chartmann et alj 
l2012l : ISaitoh et alJ2012h . 

The primary effect of the numerical diffusion associated with 
the coarse resolution in the simulations is under resolving Kelvin- 




x(R G ) 

Figure 2. Snapshot of the gas density (arbitrary units) in the r6 plane for 
the initial state of the NO model corresponding to t = 2012.5. The dashed 
yellow line shows the orbit of the cloud center-of-mass. 



Helmholtz instabilities at the cloud surface. On the other hand, the 
presence of magnetic fields in the external medium is expected to 
dampen the Kelvin-Helmholtz instability (Markevitch & Vikhlinin 
2007). Therefore, the potential effects of this instability will be less 
significant in the present problem and in a realistic disk. 

We choose an orbit orientation that allows us to utilize the 
highest possible resolution to resolve the shocks that develop in the 
cloud-disk interaction. We assume that (i) the orbital plane of the 
cloud is highly inclined (inclination angle i = 60 deg) with respect 
to the equatorial plane of the accretion flow, (ii) pericenter is in 
the equatorial plane (argumen t of periapsis oj - 0, c orresponding 
to a single disk crossing, see Sadowski et al. 2013b), and (iii) the 
cloud and the disk gas counter-rotate with respect to each other. 
This ensures that the shock normal at pericenter is nearly along 
the 6 direction, where the simulations have the highest resolution. 
As we show in the next section, the setup we use is sufficient to 
reasonably resolve the bow- shock region. For comparison, we also 
consider models with a co-rotating orbit (model N0C), two low- 
inclination models (N0LN and N0LC) and a model with oj = tt/2 
(double disk crossing, NOD). Parameters of all the models are listed 
in Tabled] 



3 RESULTS 

3.1 Physical quantities near the cloud front 

We performed 12 simulations and studied in detail the flow struc- 
ture as the cloud front sweeps past pericenter. In Fig. [3] we show 
orbital plane shapshots of density, temperature, and magnetic pres- 
sure for two representative models, NO and F0, at a time when the 
shock has reached the orbit pericenter. In both models the cloud 
orbit is inclined at an angle i = 60 deg and crosses the disk equato- 
rial plane once (oj = 0). Colors in each panel denote the magnitude 
of the relevant quantity in arbitrary units, vectors show the veloc- 
ity field, yellow dashed lines show the nominal orbit of the cloud's 
center-of-mass, and dashed vertical lines perpendicular to the cloud 
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Table 1. Cloud models 


Model 


Density 


Size 


Velocity 


^MHD 


Inclination 


Argument of 


Rotation 


Magnetic 


Relative time of 


name 


profile 




dispersion 




i 


periapsis oj 




field 


shock pericenter 


NO 


Gaussian 


FWHM=3 x 10 i:, cm 


scaled 


2012.5 


n/3 





counter- 


on 


-0.63 


Nl 


Gaussian 


FWHM=3 x 10 15 cm 


normal 


2012.5 


n/3 





counter- 


on 


-0.68 


N2 


Gaussian 


FWHM=4.5 x 10 15 cm 


no 


2012.5 


n/3 





counter- 


on 


-0.82 


FO 


Flat 


R ovA = 2.5 x 10 15 cm 


no 


2012.5 


n/3 





counter- 


on 


-0.63 


Fl 


Flat 


R out = 1.5 X 10 15 cm 


scaled 


2012.5 


n/3 





counter- 


on 


-0.65 


F2 


Flat 


R out = 3.5 x 10 15 cm 


no 


2012.5 


n/3 





counter- 


on 


-0.92 


NOF 


Gaussian 


FWHM=3 x 10 15 cm 


scaled 


2012.0 


n/3 





counter- 


on 


-0.60 


NOC 


Gaussian 


FWHM=3 x 10 15 cm 


scaled 


2012.5 


-n/3 





co- 


on 


-0.66 


NOLN 


Gaussian 


FWHM=3 x 10 15 cm 


scaled 


2012.5 


n/6 





counter- 


on 


-0.54 


NOLC 


Gaussian 


FWHM=3 x 10 15 cm 


scaled 


2012.5 


-n/6 





co- 


on 


-0.61 


NOD 


Gaussian 


FWHM=3 x 10 15 cm 


scaled 


2012.0 


n/3 


n/2 


counter- 


on 


-0.54 


NOB 


Gaussian 


FWHM=3 x 10 15 cm 


scaled 


2012.5 


n/3 





counter- 


off 


-0.70 
















-rotating 







^mhd is the epoch when the corresponding MHD simulation starts. The epoch of shock pericenter is given with respect to the center-of-mass pericenter 

crossing. 
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Figure 3. Profiles of density (left), temperature (middle) and magnetic pressure (right panel) on the cloud orbital plane at t = to - 0.63yr for model NO (top) 
and FO (bottom). Vectors show the velocity field. The yellows dashed line is the cloud center-of-mass orbit while the vertical lines show the cut used to extract 
profiles across the shock shown in Fig. [4] [6] and|7] 



front show the line along which profiles of various quantities are 
extracted and plotted in Figs.|4][6] and|7] 

The left panels in Fig. [3] show the distribution of density in 
these snapshots. The cloud has a large density and momentum rela- 
tive to the ambient gas in the accretion disk, allowing it to plow 
through the accretion disk easily. However, the tidally stretched 
front of the cloud is less massive in comparison and experiences hy- 
drodynamic and magnetic drag from the ambient gas. This interac- 
tion sculpts the front of the cloud into an oblate shape. We can dis- 



tinguish four distinct regions in the structure, which from bottom to 
top are as follows: (i) undisturbed turbulent disk gas at the bottom 
(white and brown shades, compare with Fig.[TJ, (ii) denser shocked 
disk gas above (light blue) which has been processed through the 
forward bow-shock which separates these two zones, (iii) high den- 
sity shocked cloud gas further up (black) which is separated from 
the shocked external gas by a contact discontinuity, and (iv) un- 
shocked cloud gas which lies above the reverse shock. 

The front of the bow- shock is not aligned with the orbit of 
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the cloud center-of-mass, but lies at a substantially smaller radius, 
R = 60R G instead of R = 110R G . This is because of two effects. 
First, given the initial shape and velocity distribution of the cloud, 
the parts of the cloud that reach pericenter first are those that have 
less angular momentum and hence pericenters closer to the black 
hole. Second, the interaction between the front layers of the cloud 
and the disk gas causes some additional loss of angular momentum 
in the cloud material 0. The effective area of the forward bow- shock 
is approximately A = nR 2 honV with Rf mnt ~ 40/? G - 10 15 cm This is 
consistent with the area assumed in Sadowski et al. (2013b). 

The middle panels of Fig. [3] show the distribution of tem- 
perature. The four zones previously mentioned can also be distin- 
guished in the temperature maps. As expected, the external gas is 
hot, the shocked external gas is even hotter, the shocked cloud gas 
is warm, while the unshocked cloud gas is cold. The shape of the 
bow-shock is most clearly seen in the bottom panel. 

The right two panels show the distribution of magnetic pres- 
sure. As the magnetized external gas passes through the forward 
shock, the perpendicular (to the relative velocity vector) field is am- 
plified by compression. This field then drapes around the contact 
discontinuity and contributes substantially to the total pressure in 
this region of the shocked gas. This strong field likely also inhibits 
Kelvin-Helmholtz instability at the contact discontinuity interface 
(Markevitch & Vikhlinin 2007). 

Fig. [4] shows in greater detail the one-dimensional profiles of 
various quantities of interest along the vertical dashed lines in Fig[3] 
for models NO (top) and F0 (bottom). The three dotted lines in the 
figure show the approximate locations of the forward shock (FS), 
contact discontinuity (CD) and reverse shock (RS). As can be seen, 
the shocks are not well resolved, both because the native resolution 
of the simulation is somewhat low and because the trajectory un- 
der consideration crosses the numerical grid at an angle in the 6(f) 
plane. Nevertheless, the basic features of the jumps across FS, CD 
and RS are clearly present. Across the FS, the density (green curve) 
jumps from its unshocked value up by a factor of a few, as do the 
temperature (orange) and total pressure (dark blue). The gas ther- 
mal pressure (thin blue line) dominates the pressure balance around 
the FS, whereas a region with high magnetic pressure (dashed blue 
line) is present close to the CD. The magnetic field accumulated 
here comes from the shock encounter with a region in the accretion 
flow where the magnetic pressure exceeded the gas pressure (as in 
the pre-shock region at ~ 30 Rq in Fig |4j). The jumps across the 
FS are consistent with a strong shock with Mach number of a few. 
More precisely, from the density jump pd/p u across the forward 
shock, we obtain an estimate for the shock Mach number 



M - 



2pd/Pu 



y + 1 - (y- Y)pdlp u 



2.5-3.0 



(1) 



where y = 5/3 is the adiabatic index. Alternatively, from the jump 
PdlPu in the gas thermal pressure (which dominates the pressure 



The simulations we performed do not allow us to study the impact of 
the cloud on the BH mass accretion rate. It is expected that most of the 
additional accreted material would come from "filaments" stripped from 
the cloud surface. Yet, we likely underestimate this effect because of poor 
resolution. Furthermore, the innermost region of our simulation domain is 
not self-similar (see Section 12. 3\ and therefore after scaling it up it will 
not reproduce the real Sgr A* accretion disk at the corresponding radii. 
Thus, any estimate of the cloud mass deposited at these radii is likely to be 
unreliable. 
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Figure 4. Profiles of density, temperature, velocity, gas pressure, magnetic 
pressure, and total pressure across the shock for models NO (top) and F0 
(bottom) at t = to - 0.63yr. The cloud front is moving to the right. The 
vertical dotted lines show approximate locations of the reverse shock (RS), 
contact discontinuity (CD) and forward shock (FS). 



balance in the vicinity of the forward shock, see Fig.|4j, we derive 



(y + i)Pd/Pu + y - 1 



2y 



^2-2.5 



(2) 



Both estimates concur to suggest a moderate value for the shock 
Mach number at pericenter, M ~ 2 - 3, in agreement with the as- 
sumptions of ISadowski et al.l d2Q13bh . However, it is possible that 
we slightly underestimate the shock strength because of the limited 
resolution and the resulting numerical diffusion. 

There are somewhat larger jumps across the RS, which has a 
larger Mach number since the unshocked cloud gas is very cold. In 
particular, the density of the shocked gas (between the RS and the 
CD) is quite large. As expected, the pressure is continuous across 
the CD, and correspondingly the density and the temperature have 
inversely related jumps. 

The cyan curve shows the longitudinal velocity along the cut. 
The velocity goes from that of the external gas to that of the cloud in 
two jumps, one each at FS and RS. In addition, between the shocks 
and the CD, there is a gradient in the velocity, which we discuss 
next. 

While the profiles described above are qualitatively consistent 
with the expectations for a planar FS-CD-RS double-shock sys- 
tem, the full structure we have at hand is indeed a three-dimensional 
bow-shock. This introduces some important differences, caused by 
the fact that, in a bow- shock, gas c an flow to the side and slide 
away along the contact discontinuity dLandau & Lifshitzlll95 9). In 
order to understand the differe nces specific to this ge ometry, we 
used the numerical code K0RAL JSadowski et al.l2013ab to simulate 
a toy hydrodynamical problem in which a solid spherical object 
plows through a homogeneous hot fluid with a relative Mach num- 
ber M = 2.5. 

Fig. [5] shows the results from this hydrodynamic simulation 
after the system has reached steady state. Because we assume a 
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5 6 7 8 9 10 

radius 

Figure 5. Top: Density profile and the velocity field in a toy three- 
dimensional bow shock simulation with At = 2.5 in spherical coordinates. 
Bottom: Profiles of density, pressure, temperature and velocity along the 
impact axis. 



rigid sphere, there is no reverse shock, and the surface of the sphere 
plays the role of the contact discontinuity. The forward bow -shock, 
however, is very clearly delineated in the top panel. The bottom 
panel shows profiles of density, pressure, temperature and velocity 
along the symmetry axis. The jumps in these quantities at the FS 
are as expected from standard theory of adiabatic shocks. The main 
differences from the planar shock problem appear in the post- shock 
gas. Whereas for the planar case we expect all gas properties to be 
independent of position after the shock, in the bow- shock case we 
see that the density, temperature and pressure all increase between 
the FS and the CD, whereas the velocity drops down to zero at the 
CD. These differences are caused by gas being diverted around the 
CD, as seen from the velocity vectors in the top panel. 

A comparison of the toy model results in Fig. [5] with the G2 
simulation results in Fig. |4] indicates that similar variations are 
present in the latter. The variations in the pressure and velocity, 
in particular, are quite similar to those seen in the 3D hydro dynam- 
ical problem, suggesting that the shock geometry is responsible for 
the behavior seen in Fig. [5] 

To investigate the effects of the magnetic structure of the ac- 
cretion flow on the cloud-disk encounter, we have run the simu- 
lation NOB, without magnetic fields. As compared to our fiducial 
model NO, we find that the run NOB yields a similar structure for 
the FS, which is consistent with the fact that pressure balance at the 
FS is dominated by the gas thermal component (see Fig [4]>. How- 
ever, the RS in NOB is weaker than in NO, lacking the magnetic 
support present for NO in the vicinity of the CD. 

We also study the time evolution of the shock by looking at the 
shock structure at different snapshots in the simulation. Figs.[6]and 
[7]show cuts across the FS-CD-RS structure for models NO and F0, 
respectively, at four different times. The cloud moves from left to 
right and the vertical dotted line indicates the location of the disk 
mid-plane. The various features discussed earlier are clearly seen 
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Figure 6. Evolution of density (top), temperature (middle), and total pres- 
sure (bottom panel) across the shock with time for model NO. The vertical 
dotted line shows the location of the disk equatorial plane. 



in all the snapshots and the whole pattern moves steadily to the 
right. In addition, there is evidence for some growth in the amount 
of shocked gas with increasing time. 



3.2 Pericenter crossing times 

ISadowski et all J2013bh calculated synchrotron emission from rel- 
ativistic electrons accelerated in the bow shock of G2 and showed 
that the maximum radio synchrotron emission is expected roughly a 
month after the bow shock crosses pericenter. It is evident from the 
cloud-disk interaction models we have presented here that the bow 
shock forms ahead of the front of the cloud, and reaches pericenter 
well ahead of the pericentric passage of the cloud's center-of-mass. 
We list in the last column of TableQ]the times when the cloud front 
reaches the location corresponding to true anomaly 6 = of the 
center-of-mass orbit in each of our 12 simulations. In most models 
the crossing occurs 7 to 9 months before the cloud center-of-mass' 
epoch of pericenter (to). A larger initial cloud size (models N2 and 
F2) causes the crossing to occur two or three months earlier. On 
the other hand, clouds on orbits with low inclinations (N0LN and 
N0LC) or with double crossing (NOD) are slowed down by interac- 
tion with the disk gas compared to the fiducial model NO. In these 
models, the cloud front reaches pericenter roughly a month later. 
Model NOB with no magnetic fields shows an early crossing time 
by about a month, indicating that drag due to the magnetic field 
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Figure 7. Similar to Fig.|6]but for model F0. 

draped around the CD plays a dynamically significant role in all 
the other models. 



4 DISCUSSION 

In this paper, we performed detailed magnetohydrodynamic simu- 
lations of the interaction between the gas cloud G2 and the accre- 
tion flow around the Galactic Center black hole Sgr A*. We scaled 
down the cloud to fit it within the converged region of previously 
simulated accretion flow model solved on a fixed numerical grid. 
Despite the limited numerical resolution we were able to simu- 
late the formation of the bow- shock region and qualitatively study 
its properties. We showed that for any reasonable cloud structure 
and orbit orientation that fit the current data, the shock percen- 
ter crossing takes place approximately 7 to 9 months earlier than 
the epoch when the cloud center-of-mass reaches pericenter. De- 
pending on the orbi tal model, we predic t the peak radio emission 
in Spring 2013 forlGillessen et all (120131) or late Summer 2013 for 
|PhiferetalJd2013h orbit solutions, respectively. 

The Very Large Array (VLA) carried out observations of the 
Galactic center during October 2012, December 2012, February 
2013, and March 2013 as part of a monitoring campaign]^ The data 
taken in February and March provide only upper limits to the flux 



3 Data are publicly available at 

https : //science . nrao . edu/science/service-observing. 



from Sgr A* at low frequencies because the VLA was in the low- 
resolution D-configuration during this observation and could not 
fully disentangle the fluxes of Sgr A* from that of other neraby 
sources. Intriguingly, the flux at 14 GHz, which suffers less from 
source confusion, shows a modest flux increase, which could be 
due to the extra radio emission from the bow shock or due to in- 
trinsic variab i lity of the source. Another observation was made by 
Kassimetal. (2013) using the Giant Metrewave Radio Telescope 
(GMRT) in late January and early February 2013 and resulted in 
non-detection of sub-GHz emission. 

February and March monitoring observations did not show a 
brightening but there are recent reports of brightening at 22GHz 
and 32 GHz terunthaler et al.ll2013l : lTsuboi et alJ2013h . It remains 
to be seen whether this flux enhancement is due to the bow shock. 

We encourage continued monitoring of the Galactic Center in 
the very near future. Higher resolution observations can not only 
verify a flux enhancement but may also be able to detect the dis- 
placement of the known radio source coincident with the black hole 
from the bow shock emission ~ 4400/?g away. Finally, observations 
at higher frequencies, even with the more compact array configura- 
tions, could differentiate an intrinsic brightening. Even though we 
predict the largest flux increase around a GHz, significant enhance- 
ment can be observed up to ~ 14 GHz. 
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